Stability, effective dimensions, and interactions for bosons in deformed fields 



o 
o 

Oh! 
< 
CN 



o 

CO 

-4-* 



S3 
O 

CO 

> 

CN 

o 

CO 

o 

-I— > 



i 

-a 

o 
o 



X 



O. S0rensen, D. V. Fedorov, and A. S. Jensen 
Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark 

(Dated: February 2, 2008) 

The hyperspherical adiabatic method is used to derive stability criteria for Bose-Einstein conden- 
sates in deformed external fields. An analytical approximation is obtained. For constant volume 
the highest stability is found for spherical traps. Analytical approximations to the stability crite- 
rion with and without zero point motion are derived. Extreme geometries of the field effectively 
confine the system to dimensions lower than three. As a function of deformation we compute the 
dimension to vary continuously between one and three. We derive a dimension-dependent effective 
radial Hamiltonian and investigate one choice of an effective interaction in the deformed case. 
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I. INTRODUCTION 

Bose-Einstein condensation is routinely achieved in a 
number of laboratories Q, 0, 0, EL see further descrip- 
tions in the recent monographs [a, HI • The techniques 
involve cooling and trapping of atomic gases in exter- 
nal laser fields and magnetic fields. These traps are in 
practice of cylindrical geometry 0, 0, 13 ■ For N attrac- 
tive atoms the stability criterion is experimentally estab- 
lished to be N\a s \/b t < 0.55 [3, where a s is the scat- 
tering length and bt = \f%/ (mu>) is the relevant length 
scale of the harmonic trap of geometric average frequency 
lu = 3/L0 x ujyUJ z . A reduction from three dimensions to ef- 
fectively one or two dimensions was observed experimen- 
tally [3j in the limit when the interaction energy is small 
compared to the level spacing in the tightly-confining di- 
mension. Experiments with continuous variation of the 
trap geometry from three to either one or two effective 
dimensions Q , with a two-dimensional structure 0, Il0| , 
and an effective one-dimensional geometry [Tlj | request a 
corresponding theoretical description. 

Theoretical interpretations and the underlying analy- 
ses are frequently based on model assumptions of spher- 
ical symmetry |l2l Il3| . Confinement to lower dimen- 
sions can also be studied directly without the three- 
dimensional starting point |14|. This has been done 
with a variational calculation in Gross-Pitaevskii equa- 
tion (GPE) and more recently in the GPE with vari- 
ational dimensionality ^(|. Also effects on stability of 
deformed external fields have been investigated by use 
of the GPE formulation [H [H Extreme defor- 

mations could result in effective one-dimensional or two- 
dimensional systems which can be described by effective 
interactions of corresponding discrete lower dimensions 

in, no, nn m 

In the present article we rewrite the hyperspherical for- 
mulation from reference |23| to account for a general de- 
formation of the external field. Since two-body correla- 
tions are not yet included in the wave function, this hy- 
perspherical approach resembles a mean-field treatment. 
We investigate the stability criterion in section ITTT1 Sec- 
tion IIVI contains an approach to an effective dimension 
which depends on the deformation of the external field. 



Since the interactions are presently not included in this 
effective dimension (d) , we therefore in section [3 intro- 
duce them on top of the derived d-dimensional Hamilto- 
nian. Although the choice of interactions is not unique, 
we can with some guess obtain an alternative stability 
criterion and subsequently interpret the results in terms 
of a deformation-dependent coupling strength, which is 
finally compared with known results. 



II. HYPERSPHERICAL DESCRIPTION 

A combination of magnetic fields results in an effective 
trapping potential, which can be described as the de- 
formed harmonic oscillator potential Vt rap acting on all 
the identical particles of mass m 
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where the position of the zth particle is fi = (xi,yi,Zi) 
and the frequencies along the coordinate directions q — 
x,y,z are denoted u> q . The hyperradius p is the principal 
coordinate, which is separated into the components p x , 
p y , and p z along the different axes, i.e. 
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where = rj — r*j. In the center-of-mass system the 
remaining coordinates are given as angles collectively de- 
noted by f2. 

An application here of the method presented in refer- 
ence [23] is to assume a relative wave function as a sum of 
two-particle components. In the case of a spherical trap- 
ping field each two-body component only needs depen- 
dence on p and the two-body distance ry = \/2p sin a%j 
through an angle ay. For a deformed external field it 
also needs dependence on the angle #y between the in- 
teratomic vector ry and the axis of the external field. 
The two-body component should in the cylindrical case 
then be on the form <p(p, ay , ) , which would lead to an 
angular equation in the two variables 0:12 and $12 with 
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complicated integrals. We shall here restrict ourselves 
to a wave function which is independent of hyperangles. 
This is expected to give a fair description for dilute sys- 
tems where the large distances average out all the di- 
rectional dependence (24. This is in contradiction with 
our conclusions in ref. [25| for repulsive interactions. The 
mistake was later corrected, see the treatment of attrac- 
tive interactions in ref. [23j. 

Thus, we neglect correlations in analogy to a mean- 
field treatment, so in the dilute limit the hyperangular 
average of the relative Hamiltonian is 



(H) 



H = H x + H y + H z + V 
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where d = 1 and b 2 = h/(rnu> q ) for q = x,y,z. The 
interactions Vij are averaged over all angles f2, which for 
the zero-range interaction 4ttH a s S(rij)/m for N ^> 1 
yields 
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If we replace as p x — p y = p z = p/V3, this is iden- 
tical to the average of the zero-range interaction in the 
spherically symmetric case |23j |. 

We define the following dimensionless coordinates and 
parameters: 



= Px /_2 
X ~ bj N 

bf = b x b y b z . 
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The deformation along the different axes is then de- 
scribed by f3 and 7, and s is the effective interaction 
strength. The Schrodinger equation HF(p x , p y , p z ) — 
EF(p x , p y , pz) is rewritten with the transformation 



f(x,y,z) cx {xyz) {N 2)12 F(p x , p y , p z ) 
in order to avoid first derivatives. We then obtain 

1 d 2 1 d 2 d 2 
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where e = 2NE/(hu>). Without interaction, i.e. a s = 0, 
the ground-state solution is 

f(x, y, z) = {xyz)^^/ 2 exp[-N(x 2 + y 2 + z 2 )/4] , (13) 
which for N ^> 1 is peaked at (x, y, z) = (1, 1, 1). 
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FIG. 1: Contour plots of u{x,y, z), eq. 1121 . with x = y = 
q as a function of (q, z) for s = — 0.4/3 1//6 corresponding to 
Na s /b± = —0.4 for three deformations. The values for the 
contours change by 2, 2, and 5, respectively for a) /3 = 1 
(spherical), b) /3 — 1/16 (cigar-shaped or prolate), and c) 
f3 = 16 (pancake-shape or oblate). In parts a) and b) the 
minimum at (1,1) is indicated. In part c) the descending 
path towards (0,0) is indicated. 

The wave function in eq. (|11|) is determined by the 
properties of the effective potential u. For axial symme- 
try around the z axis the x and y directions cannot be 
distinguished, that is when 7 = and /3 = b\jb 2 z with 
b\ = b x b y . This symmetry amounts to replacing p 2 and 
p 2 by p\j2 in the equations. A convenient definition for 
this case is 2g 2 = x 2 + y 2 . Equipotential contours of u 
in the (g, z) plane for g = x — y are shown in fig. ^for 
attractive interactions. For a s < (s < 0) there is al- 
ways a divergence towards —00 when (g,z) — > (0,0), see 
eq. I|12l) . However, a stationary minimum is seen in both 
figs, (spherical symmetry) and (prolate) close to 
(g, z) = (1, 1) (indicated by dots in the figure), whereas 
this minimum has disappeared for the oblate system in 
fig. and no barrier would prevent contraction (indi- 



3 



cated by arrows in the figure). For weak attraction a 
stationary minimum is present for all deformations. 

Fig. shows cuts of the potential u(g,z) along paths 
close to the bottom of the valleys (see inset). The spher- 
ical minimum (full line) is shielded by a relatively small 
barrier from the divergence for g — > 0. The minimum for 
the prolate deformation (dashed curve) is extremely sta- 
ble although the divergence for g — * still exists. For the 
oblate deformation (dot-dashed line) the local minimum 
has vanished for this attraction strength. 
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FIG. 2: The potential u(x, y, z) for s = — 0.4/3 1//6 as a function 
of q = x = y along cuts of the (g, z) plane where z = g 1 ^^ . 
The height Mbarrier °f the local maximum at top of the barrier 
is indicated for the spherical case (3=1. The inset shows 
corresponding trajectories in the (g, z) plane, compare with 
fig. Q for the three deformations. 



III. STABILITY CRITERION 

The barrier height depends on the deformation of the 
external field, see figs. ^ and |21 Extrema {xo,yo 1 Zq) of u 
in eq. (|12|l obey the three equations obtained from 
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and cyclic permutations of x, y, and z. This can be 
used to determine the critical strength s when the local 
minimum disappears. The three eqs. I|14|) can be reduced 
to find the parameter t = 1/z 2 — Zq as a function of s 
from the equation 



-2 t + 



-t 2 + 1 



(3 + 1 



q y = /3 - 7 



(15) 
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The maximum value of the right hand side of eq. I|15|) is 
reached when t is the solution to the equation 

q x t/2 q y t/2 t/2 



Vi + &l± ^1 + ^2/4 ^TTWJI 



(17) 



This solution t — i max now from eq. (|15f) gives the max- 
imum possible value of s 2 when a local minimum of the 
effective potential u still is present. Thus we have ob- 
tained the largest value of s where stable solutions exist. 

Eqs. (|14|> - (|lTf> are also obtained for the length pa- 
rameters (b x ,by,bz) of a deformed harmonic oscilla- 
tor ground state wave function by minimizing the ex- 
pectation value of the GPE mean-field hamiltonian as 
formulated by Baym and Pcthick 15]. This relates 
the present lowest-order hypersphcrical non-correlated 
results to non-selfconsistent mean- field GPE energies. 
However, the present results can be improved by includ- 
ing hyperangles in the trial wave function. This enlarges 
the variational space in ref. |15| and the solutions are 
improved. The comparison to the self-consistent GPE 
mean-field results is less direct as the present results are 
obtained in a different space and especially the hyperan- 
gular dependence provides a very different structure for 
the trial wave function. 

The results for axial symmetry (7 = 0) are shown as 
the thin solid line in fig. [3] In these units the critical 
strength s is largest for a geometry very close to spher- 
ical. Since s = Na s /b t and 6j? = b x b y b Zl this means 
that at fixed b%, or fixed volume, the scattering length 
can assume the largest negative value for the spherical 
trap. Gammal et al. |l7j | performed a time-dependent 
study with the Gross-Pitaevskii equation (GPE) which 
resulted in the critical strengths here shown as the dot- 
ted line. This is in large regions lower than the present 
result, which might be due to our neglect of the quantum 
effects of the zero-point energy, which is included later in 
this paper. The recent value for the experimental stabil- 
ity region [j| is shown as the plus and agrees with the 
mean-field model. However, the experimental error bars 
are 10% and almost includes the present result. 

These results can be compared to an analytical "spher- 
ical" approximation where the radial motion is described 
by only p while the deformed external field remains the 
same. The effective radial potential U is then obtained 
by adding centrifugal barrier and the contributions from 
zero-range interaction and the external field. The an- 
gular average replaces each of the three components p 2 
and R 2 by p 2 /3 and R 2 /3, where R is the center-of-mass 
coordinates, i.e. 

N 1 , ,2 1 , ,2 1 , ,2 

E<Wr,)> Q - \m x+ * + Z (P 2 + NR 2 ) , (18) 
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FIG. 3: The critical strength |s| = N\a a \/bt as a function of 
the deformation (3 = b\/b^ from the potential in eq. \\2\\ (thin 
solid line), from eq. 112 1 1 (dashed line), and from a mean- field 
Gross-Pitaevskii computation by Gammal et al. (dotted 
line). The thick solid line is eq. I23H obtained by considering 
the zero-point energy. The plus is the experimentally mea- 
sured value by Claussen et al. §| with the 10% error bars. 
Regions below curves are considered stable in the separate 
treatments. The double- and triple-dashed lines indicate the 
effective cross-overs to two (2D*) and one (ID*) dimensions 
from Gorlitz et al. [JJ- 

where 3^ 4 = b~ 4 + b~ 4 + b~ 4 . The stability condition 
becomes 

^=U*(A7), fc(l,0) = ^f .0.67,(21) 
,/ 3(/3 2 -t 2 ) 4 / 3 

^^'"•"t yWW^ ' <22) 

The spherical limit corresponds to 7 = and [3 = 1 
where the barrier is present when \s\ = N\a s \/b t < 0.67. 
The result for a cylinder (only 7 = 0) is shown as the 
dashed line in fig. [21 and is noticably different from, but 
numerically almost coincides with, the "deformed" treat- 
ment, thin solid line. An extreme oblate deformation 
corresponds to the two-dimensional limit where b z <C b± 
and (3 — > 00. Here eq. I|22[) yields the critical strength 
k ~ 0.4^0.6^ n/2(3~ 1 / 3 . As seen from the contour plot 
in fig. the motion is now almost confined at z = 1. 
From x = y = g we see that u{g, g, 1) only has a local 
minimum when |s| < \J 7r/2/3 -1 / 3 , which is larger than 
the value where the z motion is not fixed. This is reason- 
able since more degrees of freedom in the model lowers 
the energy. Eq. (|15|) and Baym and Pethick [15j give the 
same value in this limit, the latter obtained with a vari- 
ational study of the GPE provided that the variational 
width in the axial direction does not change due to the 
interactions. This is identical to the criterion from study- 
ing the potential u(g, g, 1), i.e. consistent with the fixed 



value z = 1. 

Analogously, in the extreme prolate limit (one- 
dimensional) where (3 — > 0, eq. ill'^ll yields the critical 
strength k ~ O^^T^E^p 1 ^ . From eq. {TSJ we ob- 
tain in this limit k ~ 3 — 3 / 4 -v/7r/5 1 ^ 6 , which has the same 
deformation dependence, but is a factor ~ 1.7 larger than 
the value from eq. (|22|l . However, fixing x = y = 1 in 
eq. i|12[l yields no critical strength since tt(l, X,z) always 
has a global minimum. Therefore, the other degrees of 
freedom are essential in this prolate limit. 

A better stability criterion can be obtained by con- 
sidering the ground-state energy E of the boson sys- 
tem, which in the non- interacting case is Eq = Ti(uj x + 
ojy + oj z )(N — l)/2 where the center-of-mass energy is 
subtracted. The system is unstable when this energy is 
larger than the barrier height t/barrier of the hyperradial 
potential U(p); see the indication in fig- El °f the corre- 
sponding height Mbanior for the reduced potential u. With 
this condition the criterion of stability is 

where 3l± 2 = &~ 2 + 6~ 2 + &~ 2 . This is seen i n fig . El (thick 
solid line) to be below the GPE calculations [lj| • The im- 
provement is here substantial compared to when the zero- 
point energy is neglected. In particular, for the spherical 
case we get N\a s \/bt — 0.53 instead of N\a s \/bt — 0.67. 
This is within the 10% error bars of the experimental 
value. The estimate of eq. Ij23(l describes the stability 
problem better than eqs. I|21|) and ll'L'l) since it includes 
the quantum effect due to the zero-point energy. An 
improvement on this would be to solve the hyperradial 
problem with interaction effects, which would lower the 
zero-point energy. The critical value would consequently 
increase, making the present value a lower bound. 

A recent variational Monte Carlo investigation of the 
stability criterion in elongated, almost one-dimensional, 
traps yielded the stability criterion tiid&id ^ 0.35 |26|. 
where tt-id ~ N/b z is the density in one dimension and 
aiD = —b±(b±/a s — 1.0326). Eq. (|2^|l can in the one- 
dimensional limit be written as N\a s \/b± < 0.66. The 
deviation between the two results might be due to our use 
of a three-dimensional zero-range interaction in this non- 
correlated model, whereas Astrakharchik et al. j2(| used 
a one-dimensional model with a zero-range interaction 
with coupling strength proportional to 1/aiD as well as a 
full 3D correlated model with hard-sphere or finite-range 
potentials. An effective potential analogous to 5{x)/<iyd 
in the general case with intermediate deformations would 
be a rewarding goal. 

According to Gorlitz et al. [3| the interaction energy 
is smaller than the energy in the tightly-confining di- 
mension when |s| < ^32/225/3~ 5/6 for the ID limit and 
when \s\ < ^32/225/3 5 / 3 for the 2D limit. These cross- 
overs are indicated by double-dashed (two-dimensional) 
and triple-dashed (one-dimensional) lines in fig. E] Since 
the critical region in each limit is below the relevant cross- 
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over, stable and strongly deformed systems can be re- 
garded as effectively one- or two-dimensional in the sense 
of these energy relations. 



The deformation of the external field effectively 
changes the dimension d of the space where the parti- 
cles move. The field changes continuously and d could 
correspondingly vary from three to either two or one. In 
order to arrive at such a description, we aim at an ef- 
fective d-dimensional Hamiltonian analogous to eq. Q 
with only one radial variable p, a deformation-dependent 
dimension d, and an effective trap length bd, i.e. 

2™H d _ 1 d d{N -x)-id_ f? 2«# 

h 2 pd(N-D-i dp P dp + b d + h 2 ' 1 

where V represents all particle interactions in d dimen- 
sions. The requirement is that the Schrodinger equation 
HdGd = EdGd with d-dimensional eigenfunction Gd and 
eigenvalue E d is obeyed, at least on average, i.e. 

J dp p d ^- l G* d {p){H d - Ed)Gd(p) = . (25) 

The lowest free solution, that is with V — 0, is given 
by eq. (|13() . In the cylindrical case we can relate the 
rf-dimcnsional function G d to this by performing the av- 
erage with respect to the angle 9 in the parametrization 
{PXiPz) = p(sm9,cos9). With inclusion of the corre- 
sponding volume elements this leads to 

f d8 cos N - 2 6 sm 2N - s 6\F(p,6)\ 2 , (26) 
Jo 

where F(p, 9) can be obtained by rewriting eqs. (|13(l and 

en. 

The characteristic energy and length can be defined by 

E d = ^(N-l), db 2 =2b 2 ± + b 2 z , (27) 

which clearly is correct in the three limits, i.e. spherical: 
d = 3 and bd = b z = bj_, two-dimensional: d = 2 and 
b± ~3> b z , and one-dimensional: d = 1 and b z 3> b±. 

In general it is not possible to find one p-independent 
set of constants (Ed, bd, d) such that HdGd = EdGd- In- 
stead we insist on the average condition in eq. 1(25) with 
Gd and Ed from eqs. (|26() and (|27|l . The result for axial 
geometry is a second-degree equation in d with one phys- 
ically meaningful root. The results for various N values 
are shown in fig. 0] The effective dimension depends on 
N for relatively small particle numbers. When N > 20, 
the curve is essentially fixed. Furthermore, the asymp- 
totic values of both d — 1 (small (3) and d = 2 (large (3) 
are reached faster for larger N since many particles feel 
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FIG. 4: The effective dimension d obtained as a function of 
the deformation parameter f3 — b\/b 2 z . Curves for larger TV 
are very close to that for iV = 20. 



the geometric confinement stronger than few particles. 
Since these effective dimensions are obtained as average 
values over p, the system might look spherical at large 
distances and strongly deformed at small distances, on 
average resulting in the curves in fig. 0J 



V. DEFORMATION-DEPENDENT 
INTERACTIONS 

The effective dimension for the non-interacting system 
possibly changes when interactions are included. The 
steps of the previous section should in principle be re- 
peated with the interactions. However, this would be 
complicated and miss the goal which is a simple effective 
Hamiltonian with a renormalizcd interaction in lower di- 
mension, see analogies in the references 0, |2(], • 

We therefore start out with a two-body contact inter- 
action with a coupling strength which is modified due 
to the deformation. This is in line with usual renor- 
malizations due to density-dependent effects [27], ; see 
also a proposed modification due to the inclusion of two- 
body correlations in the reference [2{j. So we write a 
rf-dimensional zero-range interaction with a dimension- 
dependent coupling strength g(d) as 

V d (r ij )=g(d)SW(r ij ), g(3) = ^ , (28) 

where this "d-dimensional 8 function" is defined by 
6W(r) = for r ± and J °° dr r 4 ' 1 5^ (r) = 1. The 
distance between two particles, e.g. particle 1 and 2, is 
in hyperspherical coordinates defined by ri2 = V2psina, 
where the angle a enters the angular volume element as 

dfl a = dashi*- 1 aco&^-y- 1 a . (29) 

This is valid at least for d — 1,2,3. The effective interac- 
tion V in eq. (|24|) is for N 3> 1 then given by the average 



IV. EFFECTIVE DIMENSION 
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over all coordinates except p: 

N 2 fy /2 dn a V d {V2psma) 

h 2 2N 2 (Nd/A) d / 2 a s g{d) 
2m T(d/2) p d g(3) ' 



V 



(30) 



However, this does not yield instability for d < 2 since 
the power d in p~ d is smaller than two. 

We therefore pursue another approach. Inspired by the 
forms of eqs. (|19fl and (|30|> . we write V as 



V 



h 2 2N 2 (Nd/A)P/ 2 a d 
2m~ T{d/2) ~p~P 



«3 = a s , 



(31) 



which with 0,3 = a s coincides with the result for d — 3 
if we choose p = 3. The effective potential U d in 
the cZ-dimensional Schrodinger equation corresponding to 
eq. J2U is then 



2mU d {p) _ 2N 2 {Nd/AY' 2 a d d 2 N 2 t p 2 



T(d/2) 



Ap 2 



%■ (32) 



For p < 2 this potential always has a global minimum 
and thus no collapse is present. For p > 2 there is always 
divergence to — 00 when p — ► 0. For weak attraction, 
i.e. small \a d \, there is a local minimum. This disappears 
at larger \a d \ when 



N\a d \ bY 2 2 l +P/ 2 d{ P - 2)(P- 2 )/ 4 r(d/2) 



> 



p(p + 2)(P+ 2 )/ 4 



(33) 



The criterion in eq. Ij21|) was also obtained by estimating 
when the critical point vanished. Eq. (|21|l is valid for all 
deformations, i.e. any d. In order to be able to compare 
eqs. (12 U and (|33|) . we therefore choose p > 2 such that 
eq. always is applicable. When eqs. and l(2*T|) 
agree, the effective interaction strength a d is given by 



V^(p + 2)(P+ 2 )/ 4 /3 1 /6p/rf 



If- 2 2(P~ 1 )/ 2 r(f )5 5 / 4 (;p - 2)(p~ 2 )/ 4 4 [2 + P 



(34) 



This effective interaction strength is in fig. [S] shown as 
a function of the deformation for various choices of the 
power p. The solid line shows the result for p = 3, which 
is known to be correct for = 1 (d = 3). Similarly 
the dashed line shows the result with p = d, which does 
not work for d < 2 ((3 < 0.2). Since the effective cou- 
pling strength depends strongly on the power p, we need 
further information about how the interactions enter the 
effective potential. 

An extreme deformation might lead to effectively one- 
dimensional or two-dimensional properties. Pitaevskii 
and Stringari Q collected results for the effective cou- 
pling strength in two dimensions that yields g(2) = 
/(mb z ), whereas the result from eq. I|34|) in that 
limit is larger by the factor 5 5 / 4 /4 ~ 1.9. Even though 




0.001 



1000 



FIG. 5: The effective interaction strength ad from eq. (|34l> ob- 
tained as a function of the deformation parameter /3 = 
in the large-iV limit, i.e. the connection between deforma- 
tion and effective dimension obtained from the calculation for 
iV = 20 is used for this illustration. The vertical divergence of 
the dashed line indicates the inadequacy of the corresponding 
method when d < 2. 



the results differ by a factor close to two, the right combi- 
nation of lengths shows that we have incorporated the de- 
grees of freedom in the correct manner. This was also the 
case in the previous comparison of the stability criterion 
with the one obtained by Baym and Pethick How- 
ever, as was also mentioned by Pitaevskii and Stringari 
0, in the low-density limit in two dimensions the cou- 
pling constant becomes density-dependent, which is be- 
yond the present model where correlations are neglected. 

Since p — d for d = 1 does not agree with a mean- 
ingful interpretation of the stability criterion, a one- 
dimensional system needs a different treatment. 



VI. DISCUSSION 

In conclusion, the hyperspherical method with a non- 
correlated approach yielded stability criteria as a func- 
tion of the deformation of the external field. For con- 
stant volume the highest stability was found for spherical 
traps. Effective dimensions d continuously varying be- 
tween 1 and 3 were calculated as a function of the defor- 
mation. The system can be described by a c?-dimcnsional 
effective radial potential with a d-dimensional effective 
interaction. However, this does not have an unambi- 
gious form. Applications to restricted geometries become 
simpler, where the obtained two-dimensional coupling 
strength compares reasonably with a coupling strength 
obtained by an axial average of a three-dimensional con- 
tact interaction. For the one-dimensional case an effec- 
tive coupling strength was not obtained. 

A previous approach to a d = 1 treatment by Garnmal 
et al. [3(J shows that a three-body contact interaction is 
necessary for the GPE to produce collapse in one spatial 
dimension. In the present framework a three-body con- 
tact interaction for a constant angular wave function pro- 
duces a hyperradial potential proportional to p~ 2d which 
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for any d > 1 leads to instability if the three-body cou- 
pling strength is sufficiently negative. The dotted line in 
fig. El shows the effective coupling strength for p = 2d, 
corresponding to this three-body zero-ran ge i nteraction. 

According to Astrakharchik et al. [261 |3l| a Jastrow 
ansatz for a correlated wave function and inclusion of 
two-body interactions lead to collapse in one spatial di- 



mension. According to Faddeev calculations with the 
two-body correlated model presented in reference |23j . 
a two-body interaction and inclusion of only two-body 
correlations in one spatial dimension do not lead to col- 
lapse. It seems that at least three-body correlations or 
three-body interactions are necessary in order to achieve 
a realistic description of collapse in one dimension. 
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